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We analyse the dissociation of a photon in diffractive deep inelastic scattering in the kinematic 
regime where the diffractive mass is much bigger than the photon virtuality. We consider the 
dominant qqg component keeping track of the transverse momentum of the gluon which can be 
measured as a final-state jet. We show that the diffractive gluon-jet production cross-section is 
strongly sensitive to unitarity constraints. In particular, in a model with parton saturation, this 
cross-section is sensitive to the scale at which unitarity effects become important, the saturation 
scale. We argue that the measurement of diffractive jets at HERA in the limit of high diffractive 
mass can provide useful information on the saturation regime of QCD. 
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The understanding of diffractive interactions in electron-proton deep inelastic scattering (DIS) has been a great 
theoretical challenge since diffractive events were observed at HERA There exist many attempts to describe the 
diffractive part of the deep inelastic cross-section within perturbative QCD (for an excellent review, see Ref. Q)- One 
of the most successful approaches is based on the dipole picture of DIS 0, H which expresses the scattering of the 
photon of virtuality Q 2 through its fluctuation into a color singlet qq pair (dipole) of a transverse size r^l/Q. That 
naturally incorporates the description of both inclusive and diffractive events into a common theoretical framework 
as the same dipole scattering amplitudes enter in the formulation of the inclusive and diffractive cross-sections. 
| The dipole approach revealed that the total diffractive cross-section is much more sensitive to large-size dipoles 
7-H ■ than the inclusive one Q. More precisely, it showed that unitarity, and the way it is realized, should be important 
' ingredients of the description of diffractive cross-sections, making those ideal places to look for saturation effects at 
small- a;. The saturation parametrization of the dipole scattering amplitude was quite successful in describing both 
the inclusive and diffractive structure functions. In other studies of saturation effects in diffractive DIS, nonlinear 
evolution equations for the structure function have been derived 0, 0] , new measurements proposed pH | , and fits of 
differents sets of data performed 0, . 
' In this paper, we analyse hard diffraction when the proton stays intact after the collision and the mass Mx of the 
£p-t diffractive final state is much bigger than Q 2 . This process is called diffractive photon dissociation. We extend the 
O ' study of by keeping track of the transverse momentum of the final state partons. We propose the measurement of 
the final state configuration X + jet + gap + p in virtual photon-proton collisions. In order to connect the measured 
i-Ch jet with the final-state gluon in our calculations, the jet should form the edge of the rapidity gap. The transverse 
momentum of the jet provides a hard scale necessary for the use of perturbative QCD, making our calculations valid 
even at very low values of Q 2 . 

We express the diffractive cross-section in terms of dipole scattering amplitudes, using the results derived in |13| in 



the eikonal approximation, valid at very high j*p center-of-mass energy. We show that in the context of saturation 
theory, the transverse momentum distribution of the measured jet is resonant with the scale at which the contributions 
of large-size dipoles start to be suppressed, called the saturation scale. Using the parametrization Q of saturation 
effects, we make predictions for the kinematic domain of HERA and exhibit the potential of the diffractive jet 
production measurement for extracting the saturation scale. 

The plan of the paper is as follows. In Section II we recall the derivation of ^| for the diffractive production of a 
gluon off a qq dipole. In Section III, we derive the cross-section for the diffractive photon dissociation with production 
of a gluon jet and study its model-independent properties. In section IV, we present the saturation model that we 
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use for the calculation of the jet production cross-section. Section V displays our predictions for the HERA energy 
range, and Section VI contains conclusions. 



II. DIFFRACTIVE GLUON PRODUCTION OFF A QQ DIPOLE 




FIG. 1: Diffractive gluon production off a qq dipole. xq and xi are the transverse coordinates of the quark and antiquark while 
k and y are the transverse momentum and rapidity of the measured gluon. Black points represent emission of a gluon from the 
dipole and the vertical wavy lines correspond to a color singlet interaction with a target. The amplitude for the gluon emission 
after the interaction comes with a minus sign. 



In this section we recall the derivation of 01 of the cross section for the diffractive production of a gluon in the 
high-energy scattering of a qq dipole off an arbitrary target. We shall use the light-cone coordinates with the incoming 
dipole being a right mover and work in the light-cone gauge A + — 0. In such a case, when the dipole passes through 
the target and interacts with its gauge fields, the dominant couplings are eikonal. The partonic components of the 
dipole have frozen transverse coordinates, and the gluon fields of the target do not vary during the interaction. This 
is justified since the incident dipole propagates at nearly the speed of light and its time of propagation through the 
target is shorter than the natural time scale on which the target fields vary. The effect of the interaction with the 
target is that the components of the dipole wavefunction pick up eikonal phases. 

In Fig. 1 we present the production of a gluon of transverse momentum k and rapidity y off a quark-antiquark 
dipole with transverse coordinates xq and x\. The transverse size of the dipole |a?o — £Ci | is supposed to be small in 
order to justify the use of perturbative QCD (|:Eo — x\\ <C 1/ Aqcd)- We work in a frame in which the dipole rapidity 
is not too large so that the radiation of extra softer gluons is described by quantum evolution of the target. 

The incident hadronic state is a colorless qq dipole state which has the following decomposition on the Fock states: 

\d) = \d) + \dg) (1) 
where the bare dipole \d)o is characterised by the wavefunction 

= ^2^7^\( a ^ x o)^(oi^ x i)) (2) 
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with a and a denoting colors of the quark and antiquark, and x$ and x\ being their transverse positions. The qqg 
part of the dressed dipole \dg)a is characterised by the wavefunction 



\ d 9)o = ^= J " d 2 zdy 1 -^- 



aaXa ' 



e\.(x -z) e\.(xi-z) 



\xq~z\ 2 \xi-z\ 2 



T & a I 0, x ), (a, xi), (a, A, z, y) ) , (3) 



where (a, X,z,y) characterize gluon color, polarization, transverse coordinate and rapidity, respectively. In addition, 
e\ is the transverse component of the gluon polarization vector, and T a is a generator of SU(N C ) in the fundamental 
representation. The term in brackets in @ is the well-known wavefunction for the emission of a gluon off a qq dipole 



[4j . The two contributions correspond to emission from the quark and antiquark. The only assumption made to write 
down is that the gluon is soft, that is its longitudinal momentum fraction with respect to the incident dipole is 
small. As already mentioned, we work in the frame in which only bare or one-gluon components need to be considered 
in the wavefunction \d). Softer gluons will be included through quantum evolution of the target. 

Let us denote by \t ) the initial state of the target. The outgoing state is obtained from the incoming state \d) (8> \t ) 
by the action of the S— matrix. In the eikonal approximation, S acts on quarks and gluons as (see for example 
[ISIIl): 

S(\(a,x))®\t)) = Y,[W F (x)] aa ,\(a',x))®\t) , S (\(a, \,z,y))®\t)) = £ [W A (z)] ab \(b, A, z, y))® \t) (4) 



where phase shifts due to the interaction are described by the eikonal Wilson lines W F and Wa in the fundamental 
and adjoint representations, respectively, corresponding to the propagating quarks and gluons. They are given by the 
following path ordered exponential 



W FtA {x) = Vexp^ig s J dz+ A a _ [x, z+) T£ )A J 



(5) 



with A°L being the target gauge field, and Tp,A are generators of the color group in the fundamental and adjoint 
representation. Thus, the state l^cmt) = S emerging from the eikonal interaction, is given by 



|*o»t) = 
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with 
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This outgoing state would be the one to consider to compute inclusive cross-sections, as no restrictions on the final 
state have been imposed. To compute diffractive cross-sections, one has to project the outgoing state \^ ut) on the 
subspace of color-singlet states. We have defined diffractive processes as ones in which the target does not break up, 
therefore one also has to project the outgoing state on the subspace spanned by the target state \t). Those projections 
are described in detail in [13|. They create the rapidity gap, preventing the emissions of gluons softer than the one 
described by Let us denote the resulting state \^diff), it is given by: 
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with 



I* 
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The state 1*2) represents the first contribution pictured in Fig.l when the interaction happens after the emission 
of the gluon. The second contribution when the interaction happens before the gluon emission is part of |Wf). In 
order to see that, one has to substitute \d) = \d) — \dg) in \^f) : the term that comes with \d) is the contribution of 
elastic scattering while the term that comes with \dg)o and a minus sign represents the second contribution of Fig.l. 
One can drop the elastic part since it does not contribute to gluon production and write: 
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with 



- {t\Tr(wl(xx)T a W F (x )T b ' S jWX b (z)\t) - C F {t\Tt(w^(x 1 )W F (.xo))\t) . (13) 

From this final state, one calculates the diffractive cross-section for the production of a gluon of transverse momentum 
k and rapidity y using the following formula: 

N* — l 

ijg^i) - I ^ <**//Ka(*.»)^(*.»)I**//) ( 14 ) 
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where x (k,y) and a c> \(k,y) are respectively the creation and annihilation operators of a gluon with color c, polar- 
ization A, rapidity y and transverse momentum fc. The quantity xoi =£o — x i is the transverse size of the incoming 
dipole, and b— (xq + xi)/2 is the impact parameter. The cross-section is computed in details in 0, the final result 
is: 
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with $ given by formula l|13[) . Making use of the following identity 



$(«l)**(22) 

(15) 



2Tr(^(x 1 )T a ^ F (x )T & ) [%(zf 6 = T±(wUxi)W F (z))Tr(w F (z)W F (x f) - -^Tr (w F (xi)W F (x f) , (16) 

one is able to rewrite eq. (I13fl in terms of the following S— matrices: 

S(x , Xl ) = ^-(t\T[(wUxi)W F (x ))\t) (17) 

for the scattering of a dipole with the quark and antiquark at transverse coordinates xq and x\ respectively, and 

SW(xo,z,xi) = ^(t|Tv(w4(a:i)W^W)'&(w2.(«)lVF(a;o))|t) (18) 

for the scattering of two dipoles, one with the quark and antiquark at transverse coordinates xq and z and the other 
with the quark and antiquark at transverse coordinates z and x\ . One obtains: 

^*(z) = SW(x ,z, Xl ) - ■ (19) 

We have not specified the rapidity dependence of the S— matrices, it is the rapidity at which the target is evolved. If 
Y is the total rapidity, then the S— matrices in $ depend on Y — y. 

Note that if one considers that the target is a nucleus, and that each scattering on the nucleons happens via a 
two-gluon exchange, then the target averages in l|13l) are computable (see e.g. 0) and one recovers the result of [Tlj . 
Formulae (|15|) and l|13J) are a generalization to any target that includes all numbers of gluon exchanges. Note also 
that if one writes eq. I|19[l in terms of T— matrices (S = 1 — T), one recovers the two-gluon exchange approximation 
calculated in [T^ . [l7T | by neglecting the term proportional to T 2 . Let us now apply formulae (|15(l and (|13fl to diffractive 
photon dissociation. 



III. DIFFRACTIVE PHOTON DISSOCIATION 

In deep inelastic scattering, a photon of virtuality Q 2 collides with a proton. In an appropriate frame, called the 
dipole frame, the virtual photon undergoes the hadronic interaction via a fluctuation into a dipole. The wavefunctions 
| tflp | 2 and |"0lI 2 ! describing the splitting of the photon on the dipole, are given by 

|V£(r, a; Q)\ 2 = ]T e) ((a 2 + (1 - a) 2 )e 2 Kl{e f \r\) + m 2 f K 2 (s f \r\)) 

M f ( 20 ) 

|^(r,a;Q)| 2 = «=£s £ e}4Q 2 a\l - afK 2 {e f \r\) 
71 f 



for a transversely and longitudinally polarized photon, respectively. In the above ej — y a(l — a)Q 2 +m 2 with rrif 

the mass of the quark /, r is the transverse size of the qq pair and a (resp. 1 — a) is the longitudinal momentum 
fraction of the antiquark (resp. quark). The dipole then interacts with the target proton and one has the following 
factorization 

a^ = J d 2 r j\a(\^ T (r,a;Q)\ 2 + \^ L (r,a;Q)\ 2 )a(r,a) (21) 

which relates a cross-section for an incident photon a 1 p to the corresponding cross-section with an incident dipole 
a(r,a). In the leading logarithmic approximation we are interested in, the dipole cross-sections do not depend on a 
and one defines then: 

<f>(r,Q) = J da(\^ T (r,a;Q)\ 2 + \^ L (r } a;Q)\ 2 ) . (22) 
We are going to use the factorization formula (|21|l to compute the diffractive photon dissociation cross-section. 




FIG. 2: Diffractive photon dissociation in virtual photon-proton collisions. This is the dominant contribution to the diffractive 
cross-section when the final-state diffractive mass Mx is much bigger than the photon virtuality Q. 



In diffractive deep inelastic scattering, the proton gets out of the collision intact and there is a rapidity gap between 
that proton and the final state X, see Fig. 2. If the final state diffractive mass Mx is much bigger than Q, then the 
dominant contributions to the final state come from the qqg component of the photon wavefunction or from higher 
Fock states, i.e. from the photon dissociation. By contrast, if Mx^Q, the dominant contribution comes from the qq 
component. In this paper we investigate the qqg component, in the kinematical region where 



One can easily express the diffractive mass of the qqg final state in terms of the kinematical variables of the partons: 

k 2 k- k 2 

M 2 X = -1 + + (24) 

Zq Zq Z g 

where z q , Zq, and z g are the longitudinal momentum fractions of the quark, antiquark and gluon, respectively, (z q + 
Zq + z g = l), and k q , kg, and k g are their transverse momenta (k q +kq+k g = 0). Several kinematical configurations can 
provide /3^1, however the configuration that gives the dominant contribution to the cross-section is: 

k 2 k 2 k- 

-i » A -1, Q 2 . (25) 



This is due to the infrared singularity of QCD. Indeed, as we shall see later, this ordering corresponds to the resum- 
mation of Feynman diagrams in the leading logarithmic log(l//3) approximation. It also corresponds to a final-state 
configuration where the gluon jet is the closest to the gap. 



In the previous section, we obtained the diffractive cross-section for the production of a gluon with transverse 
momentum k and rapidity y in the collision of a dipole of tranverse size xoi with the target proton in the approximation 
(J2SJ. The result reads 

-{xqi) = -7-^77- / d b A(k,x ,xi;Ai]) ■ A*(k,x ,xi;Ari) (26) 



d 2 kdy y V1J 4ir 2 C F 

where xq = 6 + xoi/2, x\ = b — XQi/2, Y is the total rapidity and At] = Y — y is the rapidity gap. The two-dimensional 
vector A is given by 

A(k,x 0) x 1 ;Ar 1 )= f e -** Uz£° _ _!z£i 1 f S ^(x ,z,x i; Ar,) - S(x 0) x i; At,)) (27) 

where S(xq, x\; Ar,) is the elastic S~ matrix for the collision of the dipole (xo, x\) on the target proton evolved at the 
rapidity At,, and (xo, z, x\\ Arf) is the elastic S— matrix for the collision of two dipoles (xq,z) and (z,x\). These 
formulae are valid at leading logarithmic accuracy in y = log(l//3) as pointed out before. Indeed, after squaring the 
term in brackets in 1)27(1. one obtains the BFKL kernel. The first (resp. second) term in the brackets corresponds to 
the emission of the gluon at transverse position z from the quark (resp. antiquark) at transverse position xo (resp. 
xi). The (resp. S) term represents the case where the interaction with the target takes place after (resp. before) 
the emission of the gluon. We shall refer to it as the real (resp. virtual) term. 

Let us introduce the usual kinematics of diffractive DIS: Y = log(l/x) and At? = log(l/xp) with 

Q 2 Q 2 + M 2 X 

W+w* ' Xr = W^' (28) 

where W 2 is the center-of-mass energy of the photon-proton collision. Using the factorization J2J, one obtains the 
qqg component of the diffractive cross-section in the virtual photon-proton collision: 

daj iff _ 2M X r ,0 ,„ , da d ^f 



d 2 kdM x mT^/V^L^M. (29) 

with the photon wavefunction given by formula (12211 and the dipole cross-section given by formulae 1|26J) and l)27|l . 
It is differential with respect to the diffractive mass Mx and to the final-state gluon transverse momentum fc, which 
can be identified with the transverse momentum of the jet which is the closest to the rapidity gap. Note that the 
transverse momentum k provides the hard scale, so that we can apply our formulae to even low values of Q 2 . 
Let us make some general comments on the dependence of the cross-section (|2l?)l . 

• When k — > 0. 

The amplitude A given by eq. (|27|l takes a constant value. The infrared divergences a priori appearing for the 
virtual term cancel between the xq and x\ part and the dominant contribution to A is determined by the large 
z behavior of (xo, z, x\\ Arf). In particular, the value of z at which {xq, z,x\] Aif) starts decreasing to 
zero plays the role of a natural cutoff and determines the value of A. The constant value of the cross-section 
(I29|l at k = is then very sensitive to the way that unitarity sets in. 

• When k — > 00. 

The amplitude A decreases as 1/fc 2 . By changes of variables, one can write 

g ik.XQ r ^2 ^ n \ Z f 1 \ \ 

A{k,x ,x 1 ;Ar ] ) = —— J — e ~ lk - z/w — (S (2) (x , z/\k\ + x , an; A77) - S(x , Xl ; A v )j - {x <-► Xl ). (30) 

Then taking k—>oo, and using 

S (2 \x , Xo ,xi; Ar/) = S {2 \x , X\, X\\ Ar/) = S(x , x%; Arj) , (31) 
one sees that the l/\k\ term vanishes leaving the dominant contribution behaving as 1/fc 2 : 



A(Mo,zi;A7j) = ^(e' tkxo V Z S^ _ e -**-*i V,£< 2 > V (32) 

K \ x xi J 

Squaring and integrating the impact parameter one obtains 

da diff . . F{\x m \) + G{\x 01 \)cos(k-x m ) 

-aWTy M ^ fc3 (33) 

with F and G depending on the precise form of S^ 2 \ When integrating over the angle of xqi in (|29|l . the G part 
becomes suppressed due to the Jo(fc|xoi|) function, and the cross-section then falls as 1/fc 4 . 



These features are general, independent of the form of the S— matrices. If one looks at the behavior of the observable 



„ dal ff 

k 'wint x < 34 > 

as a function of the gluon transverse momentum k, it is going to rise as k 2 for small values of k and fall as 1/fc 2 for 
large values of k. A maximum will occur for a value ko which is related to the inverse of the typical size for which 
the 5— matrices approach zero; in other words, the maximum kg will reflect the scale at which unitarity sets in. We 
want to explore this phenomenon in the framework of the theory of parton saturation where unitarity is realized 
perturbatively. 



IV. SATURATION MODEL FOR THE S MATRICES 



The exact form of the S— matrices is unknown, and we have to consider models in order to produce values of 
the observable (|34l) at any value of k. For this purpose we consider the following model, inspired by the GBW 
parameterization Q of parton saturation effects: 



S(x ,xi; Ar/) = 0(R P - 
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(36) 



where R n is the radius of the proton. Q s is the saturation scale, the basic quantity characterizing saturation effects 



UM LLM 1211 l2il ■ ^ i s a rising function of energy through its xp— dependence. The b— dependence of the S— matrices is 
justified if the dipole sizes \xo±\, \xq — z\ and \z— x\ \ contribute only when they are much smaller than R p . That is we 
assume that xp is always such that Q s (xp) 3> Aqcd ■ Note that the model l(3fi|) for S^ 2 ' neglects correlations between 
the two dipoles, as it is a product of two S"s. 

Interestingly enough, since the S— matrices (I35|l and H3(jfl are Gaussians, one can analytically compute the amplitude 
A given by eq. (|27(l . The details of the derivation and the final result 1A9|) are presented in Appendix A. With these 
results, we obtain for the product A ■ A*: 
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where we suppress the dependence of Q s on xp in the notation. Because of the theta function, the b— integration in 
eq. I|2t)|) gives a factor irR 2 = ctq/2. The result of the integration is then a function of Q s and the two-dimensional 
vectors k and Xqi- Inserting (|37|l into 129|l . one finally writes 
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(38) 



where now k = |fe|, and (j)(r, Q) is given by eq. 1)22(1 . The a-integration to calculate <p(r, Q) and the integrations over 
r and 9 can easily be done numerically. 



V. NUMERICAL ANALYSIS AND ITS IMPLICATIONS 



Let us analyse the k— dependence of the diffractive cross-section l|38[) . For this purpose, we define the scaled 
diffractive cross section 
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(39) 



which allows us to leave aside the problem of uncertainties due to a s and do. The inclusion of a s and ao, which have 
constant values in the kinematical domain we consider here, in the actual observable (|34|1 of course will not change 
the following discussion. In addition to the gluon transverse momentum k, cr scaled is a function of two variables: the 
photon virtuality Q 2 and the saturation scale Q s . We have to keep in mind that the diffractive cross section l|38|) was 
derived under the assumption that M\ 3>Q 2 , thus the factor in the brackets on the r.h.s. of eq. (|39[1 is close to one. 
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FIG. 3: The scaled cross section a scaled and k 2 a scaled as a function of gluon transverse momentum k for a fixed value of 
Q 2 — 1 GeV 2 and four indicated values of the saturation scale Q s - 



In Fig. 3 we plot <j scaled and J t 2 a scaled as a function of k, for fixed Q 2 = 1 GeV 2 and four values of the saturation 
scale, Q s = 0.5, 1, 2, 3 GeV. As discussed in section ITTT1 independently of the form of the S— matrix, <j scaled goes to a 
constant at small momenta while at large momenta a scaled ~ 1/fc 4 . We check that this is the case on the first plot. In 
the model with parton saturation l|35|) the value of 

^.scaled 

k — > is strongly related to the saturation scale Q s . This 
relation is better illustrated on the second plot which represents the dimensionless quantity k^a" . We see that the 
transition region between two distinct behaviours at small and large fc 2 , which features a marked bump, is linked to 
the value of Q s . It is interesting to explore this observation hoping for the possibility to extract the saturation scale 
from the measured dependence of the diffractive cross section (|34() on the gluon transverse momentum. Of course, 
in the experimental situation the gluon is seen as a jet. In the kinematic region of high diffractive mass (/3 <C 1) 
the gluon jet is the closest to the edge of the rapidity gap. The contributions from quark-initiated jets close to the 
rapidity gap in such a kinematic domain are suppressed by log(l//3). 

In order to quantify the dependence of the position of the maximum of k 2 a scaled on Q s , we plot this cross section 
as a function of the rescaled transverse momentum k/Q s for the four values of the saturation scale indicated in Fig. 3. 
In Fig. 4 we show the result of this study for two extreme values of the photon virtuality, Q 2 = 0.1 and 100 GeV 2 . 
As clearly seen, the maximum for each curve is independent of Q s and Q 2 in a broad range of considered values. 
From this figure we find that k max /Q s — 1.4—1.5, thus within the saturation model, the maximum of k 2 a scaled (k) 
is proportional to the saturation scale Q s with a coefficient of proportionality independent of Q 2 . In this way, if the 
saturation model is accurate, the diffractive gluon production in the domain of large diffractive mass offers a unique 
opportunity to determine the saturation scale Q s and its dependence on xp. 
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FIG. 4: The dimensionless cross section ] t 2 a scaled as a function of k/Q s for two extreme values of Q 2 equal to 0.1 and 100 GeV 2 . 
The four curves for each value of Q 2 correspond to four values of the saturation scale: Q s = 0.5, 1, 2, 3 GeV. 



As already discussed, in the experimental verification of the validity of our description in the ep collisions at HERA, 
one should consider large-mass diffractive processes (Mj ^S> Q) with a final-state configuration with a jet close to the 
rapidity gap: X + jet + gap + p. Then, the diffractive cross section l|34|) should be determined as a function of the jet 
tranverse momentum for different values of x-p. Positions of the maximum of the measured cross section should be 
independent of Q 2 , leading to the xp— dependence of the saturation scale. The absolute value of the saturation scale 
depends of the coefficient of proportionality between k max and Q s , which in our model equals 1.4 — 1.5. Note that 
since k max is independent of Q 2 , a wide range of photon virtuality could be used to carry out this measurament, as 
long as one keeps (3<^1. 

However, from the experimental point of view there exists an important limitation related to the minimal value 
of the transverse momentum which could be measured for a jet. In the most pesymistic scenario, considering even 
rather high values of the saturation scale, Q s (xp) ~ 1 GeV, it is unlikely that the maximum k max of the cross section 
(13411 can be seen at HERA. Thus, to see the transition between the two different behaviours of the cross section 1341) 
seems like a major experimental challenge. 

In Fig. 5 we illustrated such a situation, when the saturation scale was taken from the model Q in which 



with the following parameters: A = 0.288 (0.277) and xq = 3.04 x 10~ 4 (0.4 x 10~ 4 ), where the numbers in parenthesis 
refer to the case where, in addition to the light quarks, the charm quark is included in H20(l . The cross-section H38J1 ■ 
shown in Fig. 5, is computed for the above saturation scale and the parameters: a,, = 0.15 and Co = 23.03 (29.12) mb 
taken from [Tj- With these values and the saturation scale (|40[l . the diffractive cross-section (|34() integrated over 
transverse momentum k, da d% ^ /dMx, is well described [llj- The values of photon virtuality Q 2 , energy W and 
diffractive mass Mx indicated in Fig. 5 were taken form a recent analysis by the ZEUS Collaboration j^. One sees 
that, unfortunately, the data should always lie on the perturbative side of the bump. However, it is not necessary 
to see the whole bump to confirm the influence of the saturation scale on the results. In particular, there is a big 
difference in the rise towards the bump between the highest xp— bin (Mx =40 GeV and W — 100 GeV) and the lowest 
xp— bin (Mx = 5 GeV and W = 245 GeV). A confirmation of such a behaviour would be a sign that the saturation 
region is indeed close and could lead to the determination of the saturation scale. If however this behaviour is not 
observed, it could reflect that our saturation model is incomplete, e.g. for example Ij36|l neglects dipole correlations. 
It could also mean that in this process, unitarity does not come from saturation, but rather from soft physics |23T ]. 



Q s (x P ) = (x /xp) x/2 GeV 



(40) 



Q 2 = 0.1 GeV 2 




FIG. 5: The cross-section k' 2 Mxda / dk 2 dMx as a function of the jet transverse momentum k for Q 2 = 0.1 GeV and different 
values of diffractive mass Mx and energy W . Parameters of saturation effects are from pj: solid lines correspond to the case 
with light quarks only, while the dashed lines are predictions including the charm quark contribution. 



VI. CONCLUSIONS 



Let us summarize the main results of this paper. We analysed the contribution of the qqg component of the virtual- 
photon wavefunction to the diffractive cross sections measured at HERA. In particular, we studied in detail diffractive 
jet production in the large diffractive mass limit (/3 <C 1) with a jet close to the rapidity gap. In such a case, the jet is 
initiated by the gluon. We expressed the diffractive photon dissociation cross-section l|34|) in terms of dipole scattering 
matrices, formulae (|2bl . (|27|) and (|29[) . We found that this cross section is strongly sensitive to unitarity constraints. 
In particular, independently of the form of the scattering matrices, the cross section l|34|) is a rising (falling) function 
of the final-state gluon transverse momentum in the limit k — > (k — > oo), with the maximum related to the scale at 
which unitarity effects become important. 

In the context of saturation theory in which unitarity is realized perturbatively, the maximum is determined by 
the saturation scale Q s . Using the saturation parametrization l|35(l and l|36(l of the scattering matrices, we verify that 
the relation between the maximum of l|34JI and the saturation scale is universal, i.e. independent of Q 2 . Therefore, 
we propose the measurement of the diffractive jet production cross section in j*p collisions at HERA featuring the 
final-state configuration: X + jet + gap + p with a jet close to the rapidity gap. We argue that such a process offers 
an opportunity to extract the saturation scale from the experiment, provided a low enough jet transverse momentum 
can be measured. 
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APPENDIX A: DERIVATION OF THE AMPLITUDE A 



In this Appendix, we compute the amplitude (|27|) 

d 2 



A(k,x ,xi;x P ) = I — 



~—ik.. 



Z — Xq Z — Xi 



\z-x \ 2 \z—xi\ 



(^S (2 \xo,z,xi;x r ) - S(x Q ,xi;x r ) \ , 



with the S— matrices given by the saturation model l|35|) . The virtual contribution is proportional to 
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\z~x \ 2 \z-xi\ 2 



One can then write 

A(k,xo,xv,xp) = Q(Rp 
where we have introduced 
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Introducing 9, the angle between z and k, and cf>, the angle between r and fc, one has: 

d6 



I(k,r) = 
The angular integration gives: 
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and differentiating with iVfe, one obtains 

ik + rQ 2 s /2 
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Then performing the final integration, one gets 
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Inserting l)A8j) into (| A3|> finally gives: 
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